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We study the evolution of initially extended distributions in the coined quantum walk on the line 
by analyzing the dispersion relation of the process and its associated wave equations. This allows us, 
in particular, to devise an initially extended condition leading to a uniform probability distribution 
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Introduction. - The discrete, or coined, quantum walk 
(QW) [l| is a process originally introduced as the quan- 
tum counterpart of the classical random walk (RW). In 
both cases there is a walker and a coin: at every time 
step the coin is tossed and the walker moves depend- 
ing on the toss output. In the RW the walker moves to 
the right or to the left, while in the QW, as the walker 
and coin are quantum in nature, coherent superpositions 
right/left and head/tail happen. This feature endows 
the QW with outstanding properties, such as making the 
standard deviation of the position of an initially localized 
walker grow linearly with time t, unlike the RW in which 
this growth goes with t^/^. This has strong consequences 
in algorithmics and is one of the reasons why QWs are re- 
ceiving so much attention from the past decade. However 
the relevance of QW's is being recognized to go beyond 
this specific arena and, for example, some simple gener- 
alizations of the standard QW have shown unsuspected 
connections with phenomena such as Anderson localiza- 
tion and quantum chaos 0, |3]- Moreover, theoretical 
and experimental studies evidence that the QW finds ap- 
plications in outstanding systems, such as Bose-Einstein 
condensates [h] , atoms in optical lattices [1, 0] , trapped 
ions 8', 9j, or optical devices [lol - [l^ . just to mention a 
few. This enhances the relevance of the QW as it can con- 
stitute a means for controlling the performance of those 
systems. Apart from the discrete QW we consider here, 
continuous versions exist as well [l5|, whose relationship 
with the coined QW has been discussed in p^ . 

Surprisingly enough even the simplest version of the 
discrete QW has not been studied in all its extension. 
Specifically we refer to the fact that in almost all studies 
up to date, the initial state of the walker is assumed to be 
sharply localized at the line origin, with few exceptions. 
In [l6| it was shown that for wider initial distributions (an 
extended wavepacket), the evolution of the wavepacket is 
Gaussian-like, not showing the two characteristic outer 
peaks appearing in the probability distribution for more 
sharply localized initial conditions. In also extended 
distributions, with top-hat profile, were considered in the 
context of the superfluid-Mott insulator transition in op- 
tical lattices, but no general conclusions were drawn on 
the infiuence of these extended initial conditions on the 



long time state. It is this issue that we address in this 
Letter, and the results we obtain open the way to new 
types of distributions that the QW can exhibit, e.g. vir- 
tually flat ones, with obvious impact in applications of 
this process. 

The coined QW on the line. - In this QW the walker 
moves (at discrete time steps t G N) along a one- 
dimensional lattice of sites a; G Z, with a direction that 
depends on the state of the coin (with eigenstates R and 
L). The state of the total system at {x,t) can be ex- 
pressed in the form, 



I*:.,*) = col {R^^t,L^,t) 



(1) 



where R^^t and La;^( are wave functions on the lattice. 
As |i?a;,t|^ and |ia;,t|^ have the meaning of probability of 
finding the walker at (x, t) and the coin in state R and 
L, respectively, the probability of finding the walker at 

{x,t) is 



(2) 



and J2x^x,t — 1- The QW is ruled by a unitary map 
and a standard form is 

Rx,t+i = Rx+i,t cos 6 + L^+i^t sm 9, (3a) 
Lx,t+i = Rx-i,t sm 9 - ia;-i,t cos 9, (3b) 



where 9 G [0, 7r/2] is a parameter defining the bias of the 
coin toss {9 = J for an unbiased, or Hadamard, coin). 

The dispersion relation and the group velocity. - 
Plane wave solutions to ^ exist in the form [1^ 



, where s ~ ±, k Cz [— tt. 



exp [i {kx — w^'^^t)] 

= — arcsin (cos 9 sin k) G f i f ] 



(±) 



= N± col (cos 9 cos k ± cos t 



smf 



(4) 
(5) 



and A/± is a normalization factor (any M will have this 
meaning in the following). The dispersion relation ^ 
is represented in Fig. 1 together with the group veloc- 



ity {k) 



duj/dk associated with 
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FIG. 1: (Colour online) Dispersion relation (full line) as given 
by Eq. Q, and corresponding group velocity (dashed line). 
6 = 7r/4. 



QW group velocity has been used for determining hitting 
times [20|, and it will allow us to make simple but rele- 
vant predictions about the QW dynamics when the initial 
state is a wavepacket close to some of the eigensolutions 



$1''^ ) with f],°' a smooth 



(^) 



above, say |*x,o) 
envelope. In that case, as in any linear wave system, 
one must expect that, to the leading order, the group 
velocity govern its propagation. Hence if = ±7r/2 a 
sufficiently extended wavepacket should stay at rest be- 
cause Dg^'' (±7r/2) = 0, while if fco = it should move with 
maximum velocity Wg"^ (0) = —scosO. If the initial con- 



dition projects onto both 



'J' 1^=0/ must expect that 
the initial wavepacket splits into two, moving at opposite 
velocities given by =Fcos0. Numerical simulations of the 
QW (with Gaussian initial conditions, see below) con- 
firm these notable effects and tell us that the dispersion 
relation is a powerful tool for predicting QW dynamics 
[20j . As we demonstrate in the next sections the dis- 
persion relation @ controls not only the velocity of the 
wavepacket, but also the evolution of its shape as time 
runs, what will allow us making interesting predictions. 

Broad initial distributions: Wave equations in the con- 
tinuum limit.- The goal of this section is to find a wave 
equation for the wavepacket envelope with the help of 
discrete Fourier analysis. Given a function fx on inte- 
gers x € Z, one can define its discrete Fourier transform 
(DFT) as fk = J2x fx^~^'''^ 1 which can be inverted as 
fx = (27r)"V-J*/fee''="- Applying this DFT to the 
map ^ it is straightforward to get an explicit solution 
to the QW given an arbitrary initial condition l^'aj.o)- 



The result is I*., 



=±1 



is) 



where 



and 



2tt 



jikx-uj'-'h') 



fc,0 



, (6) 



= Ex' 



As stated, we are interested in initial conditions of 



the form \'i!xfi) = Y.s=±if'- 



envelopes fx^^ vary smoothly on x, and fco is a (car- 
rier) wave number. Then 



where the 



Z^s=±l Jk-ko 



which is peaked around fc = fco as fl!'2k„ peaked around 



fc — fco = (low frequency envelope). In such cases Eq. 
(|ni) can be written as 



=e<'^---"«^'*)F,(.,t)|<))+0(Afc), (7) 

(8) 



where Afc is the width of /j, , X = fc — fco, fl = 
and we did not modify the limits of the integral because 
of the assumed smallncss of Afc. We have introduced 
two wave functions, F± {x,t), in terms of which Px,t = 
J2s=±i \Fs {x,t)\^ + O (Afc). We let Fs {x,t) be defined 
on the reals, as there is nothing against that in Eq. ([5]), 
so that it is straightforward setting a wave equation from 
that equation. 



idtFs {x,t) = -isujidxFs - \sui2dlFs 



(9) 



after Taylor expanding VL around fco, and where w„ — 
(d"w/(ifc")j,^^Q. This equation is to be solved under the 

initial condition Fs (a;,0) — flf ^ at integer x [2^ . 

Eq. dH) is a main result of this Letter. It evidences the 
role played by the dispersion relation (jH) as anticipated: 
For distributions whose DFT is centered around some 
fco, the local variations of ui around fco determine the 
type of wave equation controlling the QW dynamics. The 
first term on the rhs gives the group velocity, already 
discussed, the second accounts for diffraction, and so on. 

Application to Gaussian initial distributions.- Two 
cases of interest of Eq. ([9]) are analyzed next, corre- 
sponding to fco = 0,7r/2 as suggested by the analysis of 
the dispersion relation. First, the case fco = yields, to 
the leading order, 

dtFs ^{s cos e)dxFs. (10) 
According to (|10p . if the (broad) initial condition projects 



onto both eigenspinors 



(±) 



two wavepackets (whose 



height will depend on the projections ^<&o^'' |^a;,o)) 
will propagate without distortion at equal but oppo- 
site velocities given by fg*"* (0) = — scos6' as commented 
above. We have checked this prediction in the original 
QW map ([3]) with initial states of the form |^'i,o) — 

AAexp — i(x/(To)^ |C), which is a Gaussian of width 
CTo- The state |C) of the coin (taken equal at any site) 
controls the projections ^^'q^'' I^'x.o)- The distortionless 

propagation at a velocity Ug'^ (0) is observed in excellent 
agreement with the prediction even for moderate values 
of CTo. Nevertheless, in all deformation of the 

wavepackets is visible after some running time, the longer 
the wider the initial distribution is. This deformation is 
controlled by the next, third order derivative term, in 
which case one has an equation similar to that derived in 
[lO| , where the role of the third order derivative was ana- 
lyzed. This type of approximation was shown to be quite 



3 



good even for localized initial conditions, where the trun- 
cation of the dispersion relation is not so well justified as 
the initial condition projects over all k values. 

The case ka — 7r/2 is more interesting for our purposes. 
It is described, to the leading order, by 



idtFs 



2tan0^^' 



(11) 



which is analogous to the Schrodinger equation as 
well as to the equation of paraxial optical diffraction, 
and the pulse propagation equation in linear optical 
fibers. The solution to under Gaussian initial 



condition is Fs{x,t) = Afexp 

where CTq is the initial width and Xq is center of the 
distribution (that remains constant as the group ve- 
locity is null in the present case), and the complex 
parameter (t) 



ity Ps{x,t) 



= 
Wcxp 



ist/aQ tanfl. 



The probabil- 



(x- 



Xof I (<TQwf 



thus re- 



mains Gaussian with 



the 



width of the distribution relative to its initial width 
CToj which grows linearly with time as soon as t > 
5(jQtan9 (note that these results are independent of s). 
Again, numerical simulations of ([3]), now with l^'^j.o) = 

A/'exp \~{x/cro) /2 + iTrx/2\ \C) are in excellent agree- 
ment with the analytical predictions. 




FIG. 2: Probability distributions for 9 — n/4 and dif- 
ferent initial conditions. In (a) the initial condition is 



-I- 



{+) 

7r/2 



fully localized at the origin, and the 



time run is t = 6 x 10^; the inset shows a magnifica- 
tion (odd sites have zero occupation probability). In (b) 

l^'i.o) = A/'exp {inx/2) sine (x/ao) |'I'^^2) is used while in (c) 

additional Gaussians (of widths ac ~ 1.1(Jo,2(jo, and Sao, 
from top to bottom) multiply the initial condition; ao = 15 
and the time run is t = 20 x 10^. 

Achieving homogeneous distributions. - In QW's a most 
desired result is that the probability distribution be as 



uniform as possible after a time. A hint towards that 
goal is given by light paraxial diffraction theory -that 
the QW follows for fco = 7r/2, Eq. ([TT]): It is a textbook 
result that the far field corresponding to a light ampli- 
tude distribution sine (x) — sin {t:x) / {t:x) is remarkably 
homogeneous within a certain spatial region (23l |. 

In order to gain insight into the problem we look 
at the solution to Eq. under an initial condition 

F {x,0) — A/" sine (x/ctq). This solution is obtainable 
by Fourier transformation of the spatial coordinate p^ . 

dkPs (k, 0) exp {ikx - ^ -^2 
-1 r+oo 



isk^t/2tant 



^ dxFs {x,0) exp {—ikx). In 
^ Y^CTo rcct (cTofc/Tr), where 



Fsix,t) 

with F, {k,0) = (27r)"^/^ 

our case Fs{k,0) — {2tt 
rect (a) = 1 if |a| < 1 and otherwise. Hence Fg {x,t) ~ 

(2^)-^^^ /+;/;;* exp [*i(/.,(fc,0], with 0,(A:,e) = 
fc^ — sk'^/2tan9 and ^ = x/t. 

As we are interested in the behavior of Fg at long 
times we can use the method of stationary phase 
[2H to evaluate the leading dependence of the inte- 
gral on {x,t). It is straightforward to get Fs{x,t) — 
w{t)~^^'^exp[-i{sTT/4 + x'^ta,n0/2t)] with w (t) = 

2tt ((To tan 9) ^ t. Hence the asymptotic analysis predicts 
a uniform probability Pg {x, t) — w (t) rect [2x/w (t)], in- 
side a segment of width w{t). Thus we should ex- 
pect, after a transient time (that can be estimated as 
t > 2(70 tan 6*), ^ distribution whose width increases 
linearly with time (its standard deviation is w [t] /\/T2). 
Inspired by this result we consider the actual QW ini- 



(+)\ 



V2/ 



tial condition l^'x.o) = A/'exp (i7rx/2) sine (x/do) 

with (To the initial width. Figure 2(b) shows results of 
the simulation of (jS]), which are in qualitative agreement 
with the discussion above: A quite uniform distribution 
is attained. For comparison, the well known result cor- 
responding to an initially localized distribution is shown 
in Fig. 2(a). Main differences are the improved degree of 
uniformity in (b) , which is free from the large outer peaks 
in (a), and the fact that in (a) even/odd sites have null 
occupation probability at odd/even times unlike in (b). 
There is however a disgusting feature in (b), namely the 
high frequency ripples that appear at the plateau. Nev- 
ertheless the situation can be improved by multiplying 
the initial condition by a Gaussian of convenient width 
(this is a kind of smoothing processing, typical of optical 
diffraction '23*1), as shown in Figure 2(c), which is another 
main result of this Letter: Almost uniform distributions 
(even reaching a top hat profile) can be obtained in the 
QW by making a judicious choice of the initial condition. 
We want to stress that these homogeneous distributions 
are so after a short transient, their homogeneity increas- 
ing with time. 

The achievement of QW homogenous distributions is a 
most desired property. From the very beginning the rela- 
tively high homogeneity of the probability distribution of 
the QW corresponding to a localized initial condition has 
been considered as a positive quality of this process for 
information purposes. In (25j the presence of some de- 



4 



coherence in the process was considered to be beneficial 
because it leads to more homogeneous distributions at a 
special time. In this sense our finding may have relevance 
as we have seen that a judicious initial condition helps 
in achieving distributions with much larger and perma- 
nent homogeneity than the initially localized case, even 
including decoherence, as Fig. 2(c) clearly demonstrates. 

Conclusions and discussion.- In this work we have 
studied the influence of the initial condition on the dis- 
crete QW on the line guided by the QW dispersion rela- 
tion and by its associated wave equations. Specifically we 
have considered the evolution of initially extended proba- 
bility distributions. We have shown that sufficiently wide 
Gaussian initial distributions propagate without distor- 
tion and small width increase at velocities that can be 
tuned, including null velocity, with a proper choice of 
the phases along the initial distribution as given by the 
phase factor exp {ikox). We have devised as well an initial 
condition that leads, after a transient, to a homogeneous 
distribution whose width increases with time remaining 
highly homogeneous at any later time. This result. Fig. 
2(c), is a main result of this paper. 



We further mention that any behavior of light diffrac- 
tion (or linear pulse propagation) can be transferred to 
the QW in the case of the resting probability distribu- 
tions (fco = its continuous limit, Eq. (|lip . 
is nothing but the paraxial diffraction equation. For 
instance, a light pattern replicates at specific planes 
when it is periodic, the so called Talbot effect [23|. In 
our case, if i^^ (a; -I- A, 0) = Fs[x,Q), where A is the 
spatial period, then |F^(x,nT)r = \Fs{x,0)\ , where 
T = (27r)" A^tane* is the Talbot period and n is in- 
teger. Other optical effects, such as those of lenses, can 
also be mimicked by introducing quadratic phase factors 
[23| in the initial condition. It is just the choice of that 
condition that can make any paraxial, linear optical phe- 
nomenon be reproduced with the QW, what can be useful 
for special implementations of this rich process. 
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